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ABSTRACT 

In this work the influence of soil mechanical properties on the displacements of mat 
foundation is studied. It was introduced the soil-structure interaction that is modeled by two 
parameters, the modulus of subgrade vertical reaction (k) and the modulus of subgrade horizontal 
reaction (27). These two parameters are dependent on the geometrical and mechanical characteristics 
of the system. It appears from this study that the modulus of vertical subgrade reaction is not an 
intrinsic characteristic but depends on the parameters of the soil and the concrete (E s v s , Eb and Vb) 
and the dimensions of the plate (so dependent on the superstructure). It is clear from this analysis 
that the foundation soil parameters are more influential than those of the plate. 

Keywords: Mat foundation, Soil- structure interaction, Mechanical properties, Finite Element 
Computation. 

1. INTRODUCTION 

Developments in civil engineering constructions and especially disorders observed in the 
supporting structures of civil engineering works pushed the designers to better take into account soil- 
structure interaction in the process of calculating the foundation structures. Thus, several authors 
have worked on the modulus of subgrade reaction that is an important parameter in consideration of 
the soil-structure interface. 

This work is a contribution to the finite element analysis of foundation slab resting on elastic 
soil. The behavioral model is established, as well as all elementary matrices of the model and the 
matrix assembly. For solving the finite element model, the numerical solver FreeFem will be used. 
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In the following, after solving the fundamental equation, the influence of different parameters 
(E s v s , Eb, Vb and e) on the behavioral model of the foundation will be highlighted. 



2. PRESENTATION OF THE CALCULATION MODEL 



For analysis of raft foundations, the soil is considered as a springs assembly (with elastic 
modulus A:) infinitely close to each other and connected by an elastic membrane (horizontal modulus 
of reaction = 2T) {Figure 1). 




7 7 7 7 7 7 ~7 2 



substratum rigide 



Figure 1 A schematisation of the problem 



The theory of plate accounting for the soil- structure interaction (biparametric model) leads to 
mat foundation behavioral law (equation 1): 

/d 4 w a 4 w d 4 w\ fd 2 w d 2 w\ . f A 

DI—-+2 - 7n , + — — - IT — — + — — + kw = q(x, y) (1) 

\dx 4 dx 2 dy 2 dy 4 ) \dx 2 dy 2 J ^ v J J v ' 



Where D is the flexural rigidity of the plate and is given by: 



E b e 3 



With: 

Eb : elastic modulus of the material constituting the plate, 

e: the thickness of the plate; 

Vb : Poisson's ratio of the plate; 

k : is the modulus of subgrade reaction. 

w : is the deflection 



The parameter k has been studied by several authors. All authors coming after Biot (1935) 
tend to give higher values to the soil reaction modulus based on input parameters (Sail, 2015). As the 
goal is to better understand the deformations, in this research (for more security), it would be better 
to use in the calculations the reaction modulus equation proposed by Biot (1935). In this equation, 
the displacements of plate points increase if the soil modulus of reaction increases. This equation is 
expressed as follow: 
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k = 



0,65Er 12 E,B 4 



l-v z 



EI 



(3) 



Equation 3 has been improved by Vesic (1963), as follow: 

l-v z Vl-v 2 )E b l) W 

where: 

E s is the modulus of subgrade, 

v s is the Poisson's ratio of the subgrade; 

B is the width of the foundation; 

Eb is the Young modulus of the concrete foundation; 

/ is the moment of inertia of the cross section of the concrete. 

T is the horizontal elastic modulus of subgrade reaction. Vlasov (1949) proposes the following 
relation: 



T = ( H 0 2 dz (5) 

To a relatively deep layer of soil where the normal stress may vary with depth, it is possible 
to use, for the function 0 (z), the non-linear continuous variable defined by equation 5a. 0(z) is a 
function which describes the variation of the displacement w(x, y) along the z axis, such that: 

0(0) = 1 and 0(H) = 0 

Selvadurai (1979) suggests two expressions of 0 (z), respectively: 

0O) = (l-§) (6a) 

0{Z) = . [ ' Ll (6b) 

H: thickness of the soil layer (depth of the rigid substratum). 
3. STUDY OF THE PARAMETER K 



Figures 2-5 present the evolution of the modulus of subgrade reaction as a function of various 
parameters of the behavioral model. These figures show a decrease in the subgrade modulus 
reaction, highly dependent on soil parameters (E s , v s ) and the geometry of the structure. For some 
fixed parameters, the modulus of subgrade reaction varies very slightly with the mechanical 
properties of concrete foundation which suggests that changes in the displacements of the structure 
will be more dependent on mechanical parameters of the subgrade that those of the concrete 
foundation . Figure 2 shows an increase in the subgrade reaction modulus with the values of the 
foundation width. For a fixed value of B (Figure 3), the modulus of subgrade reaction is strongly 
influenced by soil parameters (E s , v s ). These figures also show that for some fixed parameters, the 
modulus of subgrade reaction varies very slightly with the mechanical properties of concrete 
foundation (Figure 4), so the changes in displacements of the structure will be more related to 
subgrade mechanical parameters than to those of the concrete foundation. 
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Figure 2 - Modulus of subgrade reaction according to the Figure 3 - Modulus of subgrade reaction according to the 
plate thickness for various values of width of the plate plate thickness for various values of elastic modulus of the 

subgrade 
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Figure 4 - Modulus of subgrade reaction according to the 
plate thickness for various values of elastic modulus of 
concrete foundation 



Figure 5 - Modulus of subgrade reaction according to the 
plate thickness for various values of v s 



4. - CALCULATION USING THE FINITE ELEMENT METHOD 



4.1 - Development of the stiffness matrix for a plate element 

The type of plate element used in this research is rectangular with dimensions 2a and 2b in 
the x and y directions, respectively, and a thickness e {Figure 6). This element is called MZC 
element because it was originally developed by Melosh (1963), Zienkiewicz and Cheung (1964) as 
shown in Figure 6. The nodal displacements corresponding to each node are: 
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Wi,-^ and (1=1,2,3,4) (7) 
The function of displacement chosen for this element is: 

w = [N]{w e ] (8) 
where: 

{w e } is the nodal displacement vector containing 12 components. 

The matrix [N] containing functions of displacement (MZC rectangle shaped), is given by: 

[JV] = [N ± N 2 NM (9) 

where for each node i: 

[Nt\ = [N w N i2 ,N i3 ] for i= 1,2,3,4 (10) 



* y,i\=y/b 




Figure 6 - Rectangle finite element MZC 



According to Weaver and Johnson (1984), the shape functions are: 

N tl =1(1 + <fo)(l + VoX2 + f 0 +ri 0 -e- ri 2 ) 

N i2 = - X -b m (l + £ 0 )(1 - 77 0 )(1 + r7 0 ) 2 

iV i 3=ia^(l-f 0 )(l + r 7o )(l + e 0 ) 2 
and 

$o=$i$ ,Vo=ViV (i= 1,2, 3,4) 

where the <fj and rji are given by the following table 1: 



(ID 
(12) 

(13) 
(14) 
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Table 1 - N odal coordinates for a finite element rect angle MZC 



i 


1 


2 


3 


4 




-1 


+1 


+1 


-1 


Vi 


-1 


-1 


+1 


+ 1 



The linear differential operator generalized {d} is defined such that: 



{dx 2 ' dy 2 ' ^ dxdy}' 



and 

{0} = [d}w 

The generalized strain-displacement matrix can be given by: 

[B] = [B 1 B 2 B 3 B 4 ] 

where 



[B t ] = {d}[ty] = 



r d 2 N il 


d 2 N i2 


d 2 N i3 I 


dx 2 


dx 2 


dx 2 


d 2 N il 


d 2 N i2 


d 2 N i3 


dy 2 


dy 2 


dy 2 


2 d2Nil 


2 a 2 w i2 




- dxdy 


dxdy 


dxdy- 



(i= 1, 2, 3, 4) 



From the generalized curvatures, generalized moments can be calculated by: 
{M} = [D]{0] = [D][B]{w e ] 

The stiffness matrix of the plate element [Kp ] is given by: 
[k;] = f ne [B] T [D][B]dA = ab £ J_ 1 1 [B] T [D][B] d^dq 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



After these matrix multiplications and integrations, the stiffness matrix can be obtained for a 
finite element of plate (Zienkiewicz, 1977): 



[k x + k 2 +k 3 + k 4 ] 



12(1-V 2 ) 

where matrix kj, k 2 , k 3 and k 4 can be calculated. 



(21) 



4.2-Calculation of equivalent nodal forces 

For a plate with a distributed load (q), the equivalent nodal loads are calculated with the 
following equations (Weaver and Johnston, 1984): 
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ft = j n NjqdA (22) 
Where 

ft = abjljl 1 Nlqdfd n (23) 
After development, the load vector for an element of plate can be obtained by: 

4.3 - Determination of matrix modeling soil-structure interaction 

It is shown that the stiffness matrix [K e ] and the vector load ff e } of a finite element of a 
Kirchhoff plate can be obtained by deriving the potential energy of the internal and external forces 
acting on the plate. The rigidity of the sub-floor should be extracted from the soil deformation 
energy. For clarity, the total strain energy in the plate element and the subgrade are expressed from 
the work of Turhan (1992): 

If \d 2 w d 2 w _ d 2 w~\ T _ [d 2 w d 2 w _ d 2 w~\ , . , 1 r r , AnT , r r at i a . 

U e = -L TT(77»ZrT D — ,— ,2— - \dA + - f [w(x,y)] T k [w(x,y)]dA + 

e 2 J n e \_dx 2 dy 2 dxdyl Ydx 2 dy 2 dxdyl 2 J n e L v L v J J 1 

T 

1 r \dw dw~\ ^ \dw dw~\ 7 . 



where/ 

fi e is the area of a plate element, and all other terms have been previously defined. 

The first part of the above equation shows the conventional stiffness matrix of the 
plate [kp], differentiation of the second integral over settings (nodal displacements) returns a 
matrix[/cf ], which represents the effect of axial stress in the soil. And the last term gives a matrix 
[/cf ] which represents the effect of shearing in the soil. Thus, the total strain energy of a plate 
element can be written in the form: 

U e = \M T {[k e P ] + [fcf] + [k e t ]){w e } (25) 

Thus, the stiffness matrix for an element of the plate-soil foundation system is: 

[k e ] = [k e P ] + [k e k ] + [kft (26) 

4.3.1- Vertical bending stiffness matrix [fc^] 

The total strain energy in the soil in the vertical direction is: 

(tfft)e = \S n \w(x,y)Y k[w{x,y-)]dA (27) 

For each column of soil under the plate element, the first stiffness matrix, [fc^], is calculated 
by minimizing the total energy (Uk) e by a ratio to each component of the displacement vector {>♦>, }: 

[Kh = (28) 

Using the dimensionless coordinates (J;, ij) and combining (27) and (28) yields: 
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[kt] u = -kab-r^—f\ f\ w 2 d$dr] (29) 

L Kn ) 2 dwidwj J -l J -l * ' v ' 

Displacements at any point coordinates (g, t\) of the plate are constituted by the same form of 
functions used in assessing the stiffness matrix of a plate element. Referring to the above equations 
and substituting the shape functions of the plate member in the expression of w, give the following 
stiffness matrix (12x12) to take account of the axial stress on the soil: 

[k%] = kab j!j\[N] T [N]d^d V (30) 

This matrix was developed by Chilton and Wekezer (1990), but all coefficients are false. 
Using their coefficients, the finite element model does not provide movement when the rigid plate is 
uniformly charged. Recognizing this problem, the authors have reassessed all the coefficients of the 
matrix [fc£], and the corrected coefficients gives exactly constant q/k on the move during the 
simulation of a Winkler (1932) model for a distributed uniformly charged plate. The matrix is 
partitioned into four square matrices (6x6) : 



[k£] = kab 



sym. k k3 



(31) 



4.3.2 - Stiffness matrix due to shear deformationsfcf 

The total strain energy in the shear effect is: 

T 

Wt)e=~L h-.-H 2t — ,— \dA (32) 

v tJe 2 J n e \_dx dy\ Ydx dy\ v ' 

By minimizing this function of the deformation energy in relation to all the components of 
the displacement vector of the plate, there are given the second matrix rigidity of the foundation 
which is expressed as: 

[kfh = (33) 

Combining equations (31) and (32) and using the natural coordinates, it can be obtained the 
following equation: 

= i^i^-V-i (Vw) 2 d^ (34) 
where the displacement function, w, is as defined above. 

Substituting the expression of w in equation (34) and performing the integrations and 
differentiations, give the second matrix rigidity of the foundation, which is still a 12x12 matrix: 

= 2ta„ 1 1 (i [|] r [|] + i g]^ g]) «*, (35) 



[k e t ] = 2tab 



sym. k t3 



(36) 
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5. - DETERMINING OF THE COEFFICIENTS OF THE OVERALL MATRIX OF THE 
SYSTEM 



Using the standard procedure in the method of finite elements for the assembly of the 
elements, the global stiffness matrix is represented in the form of a half band matrix. The global 
stiffness matrix for the whole system is symbolically represented by capital letters, such as: 



[K] = EStflfcp] + [k k ] + M) (37) 

[K] = [Kp] + [K k ] + [K t ] (38) 

where N e is the total number of finite elements of plate. 
The final equation to solve is given by: 



[K]{W] ={F} 



(39) 



where [K] is the global stiffness matrix of the system, {W} the displacement vector and {FJ is the 
vector load applied to the system. For solving the finite element model, the numerical solver 
FreeFem ++ is used. 

6. PRESENTATION AND DISCUSSION OF RESULTS 

For the study using the finite element method, it is considering a lOmxlOm plate subject to a 
uniformly distributed load. The mechanical properties of soil and concrete foundation vary as 
described in the analysis by the analytical method. Figures 7-15 give the following results of the 
analysis by the finite element method for fixed parameters and zero displacements at the edges. 




OOMJStl 

am2ta 

HtHJBt 

ooxtu 




Figure 7 - Mesh of the domain by finite element 



Figure 8 - Viewing of displacements 
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Figure 9 - Isovalues displacements and deformations in 3D with mesh 



The following figures show the influence of the model parameters on the displacements of 
the mat foundation. It appears from the analysis by the finite elements method that the parameters of 
the soil (E s , v s ) have a significant influence on the behavior of the system (Figures 10 and 12). The 
elastic modulus of the concrete (Eb) has a very slight influence on the calculation model (Figure 11). 
The behavioral model of the system is almost insensitive to the Poisson ratio (v*) of the concrete 
(Figure 13). Table 2 and Figure 15 show the sensitivity of the calculation by finite elements. The 
results of the analysis by the finite elements method are comforting conclusions of analytical 
calculation. The finite elements reveal almost no lift at the edges of the plate. 



Table 2. Sensitivity of finite elements calculation 



NEF 


514 


1194 


2128 


3848 


5946 


W max{m) 


0.0345691 


0.0345881 


0.0345966 


0.0345966 


0.0345966 



NEF: Number of finite elements; w max : displacement of the center of the plate. 



Figure 16 gives the results obtained from the investigation of the influence of the ratio LIB on the 
movements of the points of the plate. These results show that the higher the ratio LIB, the greater 
w max is important. However, w max tends to a limited value for large values of LIB. 
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Figure 10 - Displacement profile following the median of 
the plate for different values of Es 



Figure 11 - Displacement profile following the median of Figure 12 -Displacement profile following the median of 
the plate for different values of E b the plate for different values of v s 
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Figure 13- Displacement profile following the median of the plate for different values of Vb 



Figure 14- Displacement profile following the median of the plate for different 

values of eb 
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Fig. 15 - Displacements of the center of the plate for 
various values of the number of finite elements 
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Figure. 16- w max for various value of L/B 



7. CONCLUSION 



After studying the influence settings on the behavioral model of the structure, the use of 
results leads to strong conclusions and remarks. Following this analysis, the elastic modulus of the 
subgrade has a very significant influence on the displacements of the concrete and the parameters of 
the concrete have almost no influence on the movement of the plate. The study also reveals firstly 
that knowledge of the mechanical properties of soil is a prerequisite for the mastery of the behavior 
of foundation structures. Also, the development of the foundation soil behavior model must 
necessarily go through a complete geotechnical characterization of materials. All results (numerical) 
regarding displacements attest that the deformations are not only function of the loads but also 
function of the mechanical properties of concrete and subgrade, which logically leads to the need to 
favor a complete geotechnical characterization of materials. 
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